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SUMMARY 

An analytical method is given for calculating the maximum wall 
temperature of a tail-pipe "burner or a ram- jet combustion chamber 
cooled by air or gas flowing through a surrounding annular passage. 
Data from step-by-step integrations of the heat flow from the com- 
bustion gas to the cooling gas, which take into consideration the 
longitudinal temperature distribution cf the combustion gases, the 
variation of heat-transfer coefficients along the burner length, 
and radiation from combustion gases to the combustion-chamber wall, 
are used to develop an empirical equation expressing the outlet 
temperature of the cooling, gas In terms of known conditions of flow 
and geometry. When the outlet cooling-gas temperature is known, 
two parameters can be evaluated that make possible the determina- 
tion of the maximum temperature of the inside wall from a working 
chart, which is a solution of the heat balance across the inside 
wall at the outlet. 

Extensive combinations and ranges of the independent variables 
(mass velocity of combustion gas, mass velocity of cooling gas, 
inlet cooling-air temperature, final combustion-gas temperature, 
burner diameter, cooling-passage height, burner length, and fins) 
have been generalized by the empirical equation, which greatly 
reduces the labor of calculating the inside-wall temperature. 

Temperatures calculated by this method show satisfactory 
agreement with corresponding temperatures measured on several exper- 
imental tail-pipe burners. Corresponding temperatures calculated 
at a station 0.7 of the distance from the flame holder to the nozzle 
outlet are within ±15° of the measured cooling-gas temperature and 
within approximately ±25° of the measured inside- wall temperature. 
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The difference between the measured and calculated inside-wall tem- 
peratures ranges from about 0.06 to 0.16 of the difference between 
the measured temperatures of the inside wall and the cooling gas. 


INTRODUCTION 

Installations of tail-pipe burners for augmenting thrust at 
take-off, in combat emergencies, and for acceleration through the 
transonic region urgently demand the solution of cooling problems 
inherent in such installations. The combustion temperatures in 
tail-pipe burners and in ram Jets reach 3000° to 4000° E. The 
durability of high- temperature alloys currently used in construc- 
tion of tail-pipe burners and ram- Jet combustion chambers is 
impaired by operating temperatures greater than about 2000° R 
(reference 1). For continuous operation, some means must be pro- 
vided to prevent overheating of the burner wall and to protect the 
aircraft structure and accessories from high temperatures. 

Some suggested methods for maintaining burner-wall tempera- 
tures within safe limits make use of stratification of combustion 
to provide cool air or gases between the wall and the hot gases 
(reference 2), a liquid- cooling Jacket (reference 3), a regenera- 
tive fuel preheater Jacket (reference 4), a ceramic combustion- 
chamber lining (reference 5), a thin wail cooled by convection and 
radiation to the free stream (reference 6), and air flowing through 
an annular cooling passage surrounding the combustion chamber 
(references 2 and 7). Cooling of exhaust pipes and tail pipes in 
which there is no combustion Is discussed in reference 7. 

Several of these cooling systems are being experimentally and 
analytically investigated at the NACA Lewie laboratory. Cooling 
by means of an annular air passage is being intensively investi- 
gated. Results of a detailed investigation on cooling of tail- 
pipe burners and ram- Jet combustion chambers are presented herein. 
The analysis is primarily concerned with heat transfer and temper- 
atures and does not consider the pressure losses through the cool- 
ing passage or the effect of the oooling system on over-all 
performance. 

Calculation of the maximum inside-wall temperature from a heat 
balance across the wall at the outlet requires the solution of- a 
fourth-degree algebraic equation. This equation Is converted into 
a nondimensional form end presented as a working chart in which 
the inside-wall temperature is given in terms of two parameters 
involving known temperatures, mass flows, geometry, and materials. 



1227 


NACA EM E9L09 


3 


The cooling-gas temperature at the outlet of the cooling pas- 
sage is not generally known, but it is a function of the heat 
received from the inside wall and the heat given up to the outside 
wall. The analysis presented herein develops an empirical method 
for obtaining the cooling-gas temperature at the outlet as an aid 
in determining the maxi mirm inside-wall temperature. Calculated 
data from a number of tedious step-by-step solutions of the heat 
balance across the inside wall are generalized in an empirical 
equation that relates the outlet cooling-gas temperature to mass 
velocity of cooling gas, mass velocity of combustion gas, inlet 
cooling-gas temperature, outlet combustion-gas temperature, burner 
diameter, cooling-passage height, and burner length. The approx- 
imate range of variables investigated in the development of the 
empirical equation are summarized in the following table: 


Burner diameter, in. 10-40 

Cooling-passage height, in 1/8-7/8 

Burner length, ft 4-10 

Mass velocity cxf combustion gas, lb/(sec)(sq. ft) . . * . 3.9-32 

Mass velocity of cooling gas, lb/(sec)(sq. ft) 1.5-51 

tfess ratio, cooling gas to combustion gas 0.015-0.33 

Total pressure upstream of flame holder, 

lb/sq. ft absolute . 550-6000 

Inlet total temperature of cooling gas, °R 400-1700 

Inlet total temperature of combustion gas, °R . . . . . 400-1700 

Outlet total temperature of combustion gas, °R .... 2000-3700 

Longitudinal fins: 

Thickness, in. „ l/l6 

Spacing, in. l/2 

Height, in l/4 and l/2 


An example problem included in appendix A Illustrates the use 
of the empirical equation and the determination of the maximum 
inside-wall temperature from the chart. 


CONFIGURATION AH3) ASSUMPTIONS 
Configuration 

An annular cooling passage formed by a concentric liner inside 
a tail-pipe burner or a ram- Jet combustion chamber Is shown in fig- 
ure 1(a). Part of the diffuser mass flow is diverted through the 
annular passage for cooling purposes. Unless otherwise stated, the 
succeeding discussion refers to a tail-pipe burner with part of the 



4 


NACA EM E9L09 


turbine-outlet gae used, for cooling, but tbe method of analysis 
applies equally well to a ram Jet cooled with air or to a tail-pipe 
burner cooled with air flowing through an external shroud. 


Basic Assumptions 

The following basic assumptions were made in calculating the 
data used to develop the empirical equation relating the outlet 
cooling-gas temperature to the independent design and operating 
variables: 

1. Steady-state conditions 

2. Flow passages of constant cross-sectional area 

3. One-dimensional flow 

4. Linear combustion- temperature rise from inlet to outlet 

5. Fully developed turbulent flow 

6. Twice nonluminous radiation from combustion gas to inside 

wall 

7. Negligible temperature drop across inside wall and no 

longitudinal conduction 

8. No heat flow through outside wall 


Discussion of Assumptions 

One-dimensional flow, - The assumption of one-dimensional flow 
is permissible in the cooling passage where tra ns verse temperature 
gradients are not large. ■ Large transverse temperature gradients 
are created in the combustion gas by the combustion pattern, which 
is largely influenced by the flame-holder design and radial clear- 
ance with the inside wall and by the manner of fuel injection. A 
uniform transverse temperature distribution is assumed in the 
absence of experimental data relating the combustion-gas temperature 
at the wall to the integrated mean combustion- gas temperature. 

This assumption is conservative because, in the usual case, the gas 
temperature at the wall is lower than the mean gas temperature; 
however, if flames should impinge directly on the inside wall, this 
assumption should still be valid. 
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Linear combustion-temperature rise. - Typical longitudinal 
distribution of the mean combustion-gas temperature is shown in 
figure 2. The data were calculated from measured static-pressure 
distributions on two 20-inch-diameter ram- jet combustion chambers 
of different length. A small simplification of this analysis is 
achieved by assuming a linear rise in combustion temperature with 
length. Trial calculations using step-by-step integration to check 
the effect of several assumed distributions of combustion-gas tem- 
perature (fig. 3) show some effect .on the outlet cooling-gas tem- 
perature and an even smaller effect on the maximum inside-wall 
temperature. 

fully developed turbulent flow. - The equivalent of fully 
developed turbulent flow is considered to exist at the inlet of 
the cooling passage and. at the inlet of the combustion chamber 
although actual flow conditions may be different. 

Badlation from combustion gas to inside wall. - Badiation 
from the combustion gas to the inside wall is difficult to esti- 
mate because of the unknown degree of luminosity of the combustion 
flames and because of inadequate data on the emissivity and 
absorptivity of luminous flames. Reference 8 states that radia- 
tion from soot in industrial furnaces is of a greater order of 
magnitude than nonluminous radiation. Luminous acetylene flames 
are reported (reference 9) to radiate roughly four times as much as 
when nonluminous. The high combustion efficiency of a burner 
results in comparatively fewer luminous particles than in a sooty 
flame and the radiant heat from the combustion gas is therefore 
probably between the nonluminous value and four times the non- 
luminous value. The calculated data used in developing the empir- 
ical equation in this report are based on the estimate of twice 
the nonluminous radiation as given in reference 8. The mass ratios 
of carbon dioxide and of water vapor to combustion-gas products 
were assumed to be the same as those corresponding to complete 
combustion cf a stoichiometric mixture of air and (same 

molecular weight as kerosene) . 

Negligible temperature drop and no longitudinal conduction. - 
The inside wall is assumed thin enough that the temperature drop 
across the wall is negligible, and the thermal conductivity and 
longitudinal temperature gradients in the wall are small enough to 
eliminate consideration of longit udina l flow of heat. 

No heat flow through outside wall. - The assumption of no 
heat flow through the outside wall eliminates an additional heat 


6 


KACA EM E9L09 


balance and greatly simplifies the calculations. The outside-wall 
temperature was therefore assumed as equal to the local cooling-air 
temperature. This assumption is satisfactory in cases where the 
outside wall is Insulated from the surroundings and gives a con- 
servative estimate of cooling requirements in cases where there is 
external heat loss. When the heat flow through the outside wall 
is appreciable, the method used for calculating the outlet cooling- 
gas temperature without heat losses can be applied and then cor- 
rected for heat losses through the outside wall. 


SIMBOLS 


The following symbols are used in this report: 
A area normal to flow direction, (sq ft) 


a,d, J,m, 
u,v,z 

B 

b 

C 

Op 

*h 

Df 

Do 

F A,E 

& 

h 

E 

k 


exponents 

W o c P,c 

UkDi 

D-Df 

cooling- passage height, — g — , (ft) 
exposed perimeter, (ft) 

specif io heat at constant pressure, (Bbu/(lb) (®R)) 
hydraulic diameter, (ft) 

outside diameter of Inside wall or burner diameter, (ft) 
Inside diameter of outside wall, (ft) 

modulus that modifies equation for radiation between black 
bodies to account for emissivities and relative geometries 

mass velocity, (lb/(sec)(sq ft)) 

heat- transfer coefficient, (Btu/(br)(sq ft)(°R)) 

factor in equation (19) 

thermal conductivity, (Btu/(hr) (sq ft) (°E/ft) ) 
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L 

l 

n 

P 

P 

a 

s 

s 

T 

t 

U 

T 0 

w 

x 

y 

a 


P 

€ 

€ *W 

a 


"beam length, (ft) (reference 8, p. 69) 
length or height, (ft) 
number of fins 

total pressure, (lb/sq. ft absolute) 
static pressure, (lb/sq. ft absolute) 
heat flow, (Btu/hr) 
surface area, (sq. ft) 

distance between surfaces of adjacent fins, (ft) 
total temperature, (°E) 

arithmetic average of surface temperature and free-air 
temperature, (°E) 

static or surface temperature, (°R) 

effective over-all heat-transf er coefficient 

velooity of undisturbed air stream, (ft /sec) 

mass flow, (Ib/sec) 

distance, (ft) 

thickness, (ft) 

absorptivity- 

constant of proportionality- 

density of air at temperature T m and prevalent 
pressure, (lb/ou ft) 

emissivity 

pseudoemissivity of wall 

Stefan-Boltzmann constant, 0.173X10“® (Btu/(sq. ft) (hr) (°R^) ) 
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Parameters : 


cp 


^g(P A ,E + e V a g) 


h g + b c (l - s[) + ^ V h ° Ck f y f z f 


f ^oC 

Vf 


T g 4 ° £ Vg + Vs + + T c 


♦ 


hoi 1 


-3?) + Ji *** *f 


hg + h 0 (l - ^ ^ ^/hoC Vf tanh l f 


Subscripts: 
ay average 

CO g oarbon dioxide in combustion gas 

0 cooling gas or cooling passage 

d diffuser outlet 

f fin 

g combustion gas or combustion chamber 

H 2 O water vapor in combustion gas 

1 inside wall 

o outside wall 

r radiation 

w wall 

x variable station 

0 ambient or surrounding 

1 inlet station or preceding station (fig. 1) 

2 outlet station or succeeding station (fig. 1) 
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HEAT BALANCE AND DETEEMENATION OF WALL TEMEEEATUEE 

Under equilibrium conditions, the heat transfer from the com- 
bustion gas to the inside wall is balanced by the heat lost to the 
coo ling air, and the equilibrium temperature of the wall can be 
obtained from the heat-balance equation in terms of known tempera- 
tures, heat- transfer coefficients, emissivities, and geometry. 
Inasmuch as the heat balance is a fourth-degree algebraic equation, 
a chart is developed that reduces the labor of calculating the 
inside-wall temperature. 


Heat-Transfer Equations 

Convective heat transfer. - Convective heat transfer between 
a gas and a surface is given by Newton’s law of cooling 

g >= b(Tg - ty.) (l) 

where h depends on the fluid flow conditions. The convective 
heat-transfer coefficient h for turbulent flow of air in long 
ducts may be obtained from reference 10 as 


h 


0.378 


giO . 3 q0 . 8 



( 2 ) 


Equation (2) is used to estimate convective heat-transfer coef- 
ficients in the combustion chamber and the cooling passage. 

For circular ducts, the hydraulic diameter is equal to 

the geometrical diameter. The hydraulic diameter of a noncircular 
duct is four times the flow area divided by the wetted perimeter. 
In the case of an annulus, the hydraulic diameter is 

Dh = D 0 - (3) 

and for a rectangular fin passage 

2s Z- 

w 

where Zj is the fin height and s is the distance between sur- 
faces of adjacent fins. 
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Badlant heat transfer from combustion gas, - The radiant heat 
transferred from the combustion gas to the Inside wall is given 
(reference 8) by the equation 

^ “ oe 'w (eg T g 4 “ a g^i 4 ) ( 5 ) 


where 


c w +1 


w 


( 6 ) 


and 


*6 


and 


a S 


are evaluated from charts in reference 8. Emis- 


sivities of surfaces are given in references 8 and 11. 


Radiant heat transfer between walls of cooling -passage. - If 
the outside-wall temperature is assumed equal to the local cooling- 
gas temperature, the radiant heat transfer between the walls of the 
cooling passage is expressed by 



<**A,E ^i 4 " ^o 4 ) 


( 7 ) 


where E. _ is a combined shape- emissivity modulus that modifies 

Ai fJSl 

the equation for radiation between black bodies to account for the 
emissivities and the relative geometries of the radiating surfaces. 
Reference 10 gives E^ -g for infinite concentric cylinders as 



( 8 ) 


which is sufficiently accurate for the annular cooling passage with 
or without longitudinal fins in the passage. 

Heat transfer to fins. - The heat transfer from a finite bar 
fin (fig. l(b}) of constant cross-sectional area with no heat flow 
through the end of the fin tip is given in reference 12 as 



q^ = hoCkj. Af (tf - T c ) tanh If 


( 9 ) 


L2ZT 


where 


Af ® y f per unit of length 


Use of equation. (9) In this analysis la limited to oases -where fins are spaced close 
enough to maintain an essentially constant wall temperature "between fin "bases. 


Heat Balance 

The balance of heat received and given up by the Inside wall Ms 

lr,g + % " <*r,v + 4c + »lff 

Substitution of equations (l) , (5 ) , (7), and (9) into equation (10) and rep- 
resenting S per unit of length gives 

^VVs 4 - Vi 4 ) <®I + V T S - *i> **■ 

' or A,s(' l i 4 ' T c 4 ) ”% + M*! “ T oH’ a) l - ”?f) + nyb o 0tfy f (tj - T„) tenh t f , 
Dividing equation (11) by itD^ and expanding and collecting like terms give 


t i 4 a ( F A,E + e V a g) + H h g + K (l “ V h o Ck f7f 1 

= T g 4 ^V € g + + T c^ r A,E + T o h o (*■ “ V*V^ 



? y f tanh ^ yjj— (12) 
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Determination of Inside -Wall Temperature 


The unknown inside-wall temperature t-j_ is related in equa- 
tion (12) to the known temperatures, mass flows, geometry, and 
materials. Rearranging equation (12) gives 


<T> 




(13) 


where 


Tg 4 ° £ Vg + Vg + t ° 4 ° f a,e + *o 


* 


h c(l - Sj) + ^V b c C ^f tanh Zf 


hg + h 0 (l - ^ J + ■— V h ° C1 ^ y f ^ 




(14) 


and 



With no fins in the cooling passage, n is zero and the parameters 
\J/ and cp reduce to 


\l/ = 


= -S. 


oc 


'„e„ + T CT h CT + T^ 4 gF a>e + Tfvh, 






L C“C 


kg + k c 


(16) 


cp = 


\J/ 3 a(F 


A^jE 


+ £' a ) 
w g' 


h g + h c 


(17) 


The temperature of the inside wall can now be read from a 
plot (fig. 4) of equation (13) when the parameters cp and ^ are 
evaluated. It is noted that cp contains Og, which depends on 
the value of t^ and must therefore be estimated. At a wall tem- 
perature of 1960° E, the average rate of decrease of otg with t^ 
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Is about 10 percent per 100° increase in t^ and, because the con- 
tribution of e’v^ to od is small, a second estimate of ctg Is 
generally unnecessary. She outlet temperature of the cooling gas 
that appears in the parameter '4' is obtained from an empirical 
equation developed in the following section. 


EMPIRICAL EQUATION FOR COOLING-GAS TEMPERATURE 

Direct solution of the differential equation relating the out- 
let cooling-gas temperature to the independent variables, passage 
height b, passage length Z, burner diameter D-j., mass velocity 
of combustion gas Gg, mass velocity of cooling gas G c , 
combustion-gas temperature distribution, and inlet cooling-gas 
temperature, is impossible because the heat-transfer coefficients 
vary with temperature along the length of the burner. The dif- 
ferential equation for the temperatures in the cooling passage can 
be solved when the simplifying assumptions of constant heat-transfer 
coefficients and constant gas properties are made. Solution of the 
differential equation using these simplifications gives the charac- 
teristic equation 


T g,2 ~ t c, 2 
T g,2 " T g,l 


1 - 

P T gjl " 

/ zu N] 

ZU 

| bG c Cp^c 


- T S.3- T g.l 

\ bG c°p,cj 



ZU 


^ G c c p,c 


(18) 


which Is derived in appendix B and plotted in figure 5 for given 


T 


values of the dimensionless parameters 




-T 




T g,2‘ T g,l 


and 


ZU 


bG c 


c p,c 


Fitting Empirical Equation 

i 

The succeeding paragraphs develop an empirical equation, which 

evaluates the parameter — from the results of a number of 

tG c c p,c 

step-by-step solutions for the outlet cooling-gas temperature. The 
data obtained from the step-by-step solutions are presented in 
table I. 
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t> c 

Hi© unknown quantity in the parameter is the effec- 

tive over-all heat-transfer coefficient U. The value of U is 
related to Reynolds number of the cooling-gas flow, the Reynolds 
number of the combustion-gas flow, and the temperature of cooling 
gas and combustion gas. Because of this relation, the parameter 
bG c c c 

— Uj 2, ■ is related to the independent variables by the following 
approximate empirical equation: 


bG c c. 


Ri 


ZU 


= K 



Gg^Jb^Z^Tg^/Tc,! 2 


(19) 


The constants in this equation can be adjusted to fit the results 
of the step-by-step solutions that take into consideration the 
variation in heat-transfer coefficient. This fitting is accom- 
plished by selecting the proper values for the factor K and the 
seven exponents a, d, J, m, u, v, and z. 


The method of selecting the values of the exponents is illus- 
trated by the selection of a. From the results of the step-by- 
step calculations (table I), in which G c /Gg was varied and all 


other quantities held constant, values of 


- o,2 


The corresponding values of 






T g,2 _T g,l 


are obtained. 


T g,2 _T g,l 


are found from the initial 


conditions. The values of 


bGr.0. 


ZU 


obtained from figure 5 are 


then plotted against G c /G g on logarithmic coordinates, as shown 
in figure 6(a). The data fall in straight lines and each line 
represents a constant set of operating conditions and geometry. 
The slope of each line is the exponent a. The value of a in 
equation (19) is assumed to be independent of the variables,’ this 
assumption is shown to be approximately correct by the fact that 
lines representing widely different conditions are practically 
parallel. 


A similar procedure is employed in determining the values of 
the other exponents (figs. 6(b) to 6(g)). The extent to which 
the value of each exponent is independent of the values of the 
other variables can be seen by the approach to parallelism in fig- 
ures 6(b) to 6(g). The values of the exponents thus determined 
are as follows: 


LZZl 
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a = 0.84 
d = 0.29 

3 = 0 

m = 1.2 
u = -1.0 
v = -0.38 
z = -0.23 


The values of the variables 
and the corresponding values of 

stituted into equation (19) along -with the tabulated values for 
the exponents, to determine the factor K. The mean value of K 
(fig. 7) is 50,000 with a spread of ±0.14 from the mean value, 
which is an indication of the over-all aocuracy of the empirical 
equation in correlating the range of combinations investigated in 
the step-by-step calculations. 


from the step-by-step solutions 
from figure 5 are sub- 


^c°p , c 


In a ram Jet having an outlet combustion-gas temperature of 
3500° E, a combustion- chamber-inlet temperature of 500° E, and an 
inlet cooling-air temperature of 500° E, an uncertainty of ±0.14 
in K causes variations in cooling-air temperature rise frcm the 
mean temperature rise of the cooling air of about ±0.15 and ±0.03 
W C°B t O 


at 


ZU 


equal to 1.5 and 5, respectively. In these regions 


, the effect of outlet cooling-air temperature on the maximum inside- 
wall temperature is approximately 1° change in wall temperature per 

bCr-c,, c 

10° change in cooling-air temperature for — - TT ^ 7 -- of 1.5 and 

tu bG c c p,c 

approximately a 1° per degree change in the region of — • ■ 
of 5. 


In tail-pipe burners cooled by air entering a cooling passage 
at 500° E and having combustion-gas temperatures of 1700° and 

r\ bGrfjC'Q r* 

3500° R at the inlet and outlet, respectively, — Tn~ ■* 3as a 

value of approximately 5, and the possible variation in cooling- 
air temperature rise from the mean value of cooling-air temperature 
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rise is approximately ±0.15. For this region, a 4° change in 
outlet cooling-air temperature results in a 1° change in the maxi- 
mum inside-wall temperature. 

Rewritten in a working form, equation (19) is 


bG d 0 p,c 


ZTJ 



, 0.29 b 1,2 / 1 \ 0,38 



0.23 

) 

( 20 ) 


Figure 5 in combination with the desired empirical equation (20) 
relates the outlet cooling-gas temperature to the independent 
variables . 


Effect of Longitudinal Fins in Cooling Passage 


The applicability of equation (19) has been investigated for a 
few combinations having fins attached to the burner wall. The 
exponents a and m in equation (19) were determined and found 
to be in agreement with the exponents in equation (20). The fac- 
tor K in equation (19) was then determined using the remaining 
exponents of equation (20) . 


The form of equation (19), as indicated by figure 8, is 
applicable to longitudinal fins in the cooling passage when K is 
replaced by a function of the ratio of heat- transfer surface on 
the cooling -gas side including fins S Q to the heat-transfer sur- 
face on the combustion-gas side Sg of the Inside wall. Figure 8 
yields the approximate expression 


K =6100 




( 21 ) 


when the cooling passage contains fins of l/l6 inch (0.0052 ft) 
thick Inconel having heights of l/4 to l/2 inch (0.02083 to 
0.04167 ft) and spaced l/2 inch apart. A more thorough analysis 
is recommended for other combinations of materials and fin geometry. 


Correction for External Heat Losses 

Flight or test installations of tail-pipe burners that are 
not insulated may have appreciable heat loss through the outside 
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vail. In these cases, an approximate correction can he applied to 
the outlet cooling-gas temperature (obtained from equations (20) 
end (18)) before substitution into the parameter \J/. The heat 
loss from the outside wall to surroundings is 

= c3?A,E^o,av “ ^0 5 + ^0^0,87 “ ^O^ (22) 

where t Q aT Is the average outside-wall temperature and h D Ib 
the outside-wall heat-transf er coefficient, which may be for free 
convection or forced convection depending on the installation. 

The heat-transfer coefficient for free convection about a horizon- 
tal cylinder Is given in reference 8 as 


^o 



(23) 


If the engine and tail-pipe burner are in an unconfined air stream 
of high velocity, ho can be estimated from the flat-plate equa- 
tion (reference 10) for turbulent flow. 


h Q = 0.64 T, 


0.3 


(VqP) 


0.8 


•m 


, 0.2 


(24) 


where l is the length of the plate. The reduction in outlet tem- 
perature of the cooling gas due to heat losses to the surroundings 
Is 


AT C 


J 3600 W 0 c p ^ 0 


(25) 


and the corrected outlet cooling-gas temperature equals the outlet 
cooling-gas temperature without losses minus AT C . 


COMPARISON OF ANALYSIS WITH EXPERIMENT 

Temperatures of the cooling-passage walls and the cooling gas 
have been measured on several tail-pipe burners during an inves- 
tigation of their performance and operational characteristics in 
the NACA Lewis altitude wind tunnel. The method by which these 
data were obtained did not permit a direct check of the values of 
the factor K and the exponents in the empirical equation. A 
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preliminary check has been made, however, hy comparing the exper- 
imentally observed temperatures with the corresponding temperatures 
calculated using the empirical equation. Typical results are pre- 
sented for one of the tail-pipe burners, which is shown in fig- 
ure 9, cooled by part of the turbine-outlet gas flowing through 
the annular cooling passage and by heat losses through the unin- 
sulated outside wall (figs. 10 and 11). The heat losses to the 
tunnel test section were great enough that the outlet cooling-gas 
temperature was approximately equal to the inlet temperature. 

The outlet temperature of the cooling gas obtained from the empir- 
ical equation is therefore corrected for heat losses to the tunnel 
by the previously described method. In the calculations, the tun- 
nel is considered to be a black body, the emissivity of the inBide 
and outside walls is 0.7, nonluminous radiation is assumed from 
the combustion gas to the inside wall, and the characteristic 
length in equation (24) is 20 feet. (A sample calculation is given 
in appendix A.) The corrected oooling-gas temperature is used in 
determining the inside-wall temperature. 

Comparison is made of the calculated and measured inside-wall 
and cooling-gas temperatures at the outlet of the 48-inch-long 
cylindrical section of the combustion chamber because temperatures 
were not measured in the nozzle section of the cooling passage. 
Inasmuch as the combustion-gas temperature at the outlet of the 
48-inch section was not measured and because combustion continues 
in the nozzle, the combustion-gas temperature was estimated to 
equal the inlet combustion-gas temperature plus 0.85 of the 
combustion-temperature rise to the nozzle outlet. 

Besults and discussion. - Calculated and measured temperatures 
of the cooling gas and the inside wall are shown in figures 10 and 
11. The corrected cooling-gas temperature is within ±15° of the 
measured cooling-gas temperature, and the calculated inside-wall 
temperature is within approximately ±25° of the measured inside- 
wall temperature. The difference between the measured and calcu- 
lated inside-wall temperatures ranges from approximately 0.06 to 
0.16 of the difference between the measured temperatures of the 
inside wall and the cooling gas. 

The good agreement shown is due in part to the choice of 
combustion-gas temperature at the 48-inch station and the emis- 
sivity of the walls and of the ratio of luminous radiation to non- 
luminouB radiation from the combustion gas to the inside wall. 

The effect of changes in the chosen values upon the agreement 
between calculated and measured temperatures will be subsequently 
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discussed. The agreement between Title calculated and measured tem- 
peratures applies only to combustion chambers having a combination 
of a well-proportioned flame holder and an optimum fuel- injection 
system, which completely fills the burner cross section with flame 
and produces high combustion efficiencies. The analysis predicts 
cooling-gas and inside-wall temperatures that are too high for 
configurations in which fuel is concentrated near the center of the 
combustion chamber, or which have no flame holder, or when the 
radial clearance between the annular flame holder and the inside 
wall is large enough to provide a relatively cool layer of gas 
next to the inside wall. 

Effects of uncertainties in calculated temperatures. - The 
accuracy to which the temperatures of the cooling gas and the 
inside wall can be calculated depends on the accuracy to which the 
quantities in the empirical equation and the parameters cp and 
are known and in turn upon knowledge of actual conditions exist- 
ing in the combustion chamber and the cooling passage. 

The effects of several combinations of assumptions regarding 
combustion-gas temperature, emissivity, and ratio of luminous to 
nonluminous radiation on the computed temperatures at the outlet 
of the 48-inch section of the tail-pipe burner are shown in the 
following table. The measured quantities were: mass velocity of 

combustion gas, 6.83 pounds per second per square foot; mass veloc- 
ity of cooling gas, 9.14 pounds per second per square foot; mass 
ratio W c /W"g, 0.088; combustion-gas temperature at inlet, 1704° R; 
combustion-gas temperature at nozzle outlet, 3363° R; cooling-gas 
temperature at outlet of 48- inch section, 1712° R; average outside- 
wall temperature t Q av , 1513° R; and inside-wall temperature ^± } z> 
2130° R. The outside-wall heat- transfer coefficient h Q (equa- 
tion (24)) is estimated from turbulent flow over a flat plate 
20 feet long. 
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* Meas- 
ured 
condi - 
■ tion 

Assumed or calculated condition 

Combustion-gas tem- 
perature, T g 2> at 
48-inch station, 

°E 

i 

3363 a i 

3114 to 

Theoretical outlet 
cooling-gas tem- 
perature, °E 
(equations (20) 
and (18)) 


! 

1853 1 

1830 

Emissivity of 
cooling-passage 
wall, e w 


0.7 

0.8 

0.7 

0.8 

Corrected cooling- 
gas temperature, 

°B 

(equations (25), 
(22), and (24)) 

1712 

1734 

1724 

1711 

1701 

Eadiation from gas 
to wall 


m 

(4) 

i 

(o) 

(4) 

(°> 

(4) 

( c > 

( d > 

Inside-wall temper- 
ature, ti^gj °R 
(equation (13)) 

2130 

2215 

2315 

2188 

2290 

2118 

2195 

2094 

2185 


Combustion-gas temperature at nozzle outlet. 

^ Inlet combustion-gas temperature plus 0.85 combustion-temperature 
rise to' nozzle outlet. 

Conluminous radiation from combustion gas to inside vail. 

“Two times nonluminous radiation from combustion gas to inside 
wall. 

The preceding table shows that either of the assumed tempera- 
tures for combustion gas at the outlet of the 48-inch section gave 
corrected cooling-gas temperatures close to the measured tempera- 
ture of 1712° E. The temperatures of the cooling gas and of the 
inside wall are approximately 23° and 100° higher, respectively, 
when the temperature at the nozzle outlet is assumed at the 48-inch 
station. 

An increase In e v from 0.7 to 0.8 lowers the outlet cooling- 
gas temperature 10° E and lowers the inside-wall temperature 10° to 
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27° B. It Is therefore concluded that the choice of wall emis- 
sivity is of secondary consideration in the region of c v equal 
to about 0.7 or 0.8. Doubling the nonluminous radiation from the 
combustion gas increases the inside-wall temperature from about 
80° to 100° E. 

Obviously, a number of combinations of assumptions, within 
certain limits, will yield calculated temperatures equally close 
to the measured .temperatures. The proper choice, however, depends 
on more precise knowledge of the actual conditions. Uncertainties 
in the actual combustion-gas temperature near the inside wall are 
most critical and should be experimentally investigated for various 
combinations of flame holder and fuel- injection system. The ratio 
of luminous radiation to nonluminous radiation from the combustion 
gases is of equal importance, although more difficult to determine 
experimentally. 


SUMMARY OF BESULTS 


An analysis has been made that reduces the labor of calcu- 
lating the maximum inside-wall temperature of the combustion cham- 
ber of a tail-pipe burner or a ram- jet combustion chamber cooled 
by air or gas flowing through a surrounding annular cooling passage. 
The maximum temperature of the inside wall is obtained from a heat 
balance across the inside wall at the outlet after the outlet tem- 
perature of the cooling gas has been determined from other known 
conditions. The calculation of the inside-wall temperature from 
the heat balance requires the solution of a fourth-degree algebraic 
equation, which Is simplified by means of a working chart. 


Because the outlet temperature of the cooling gas T c> % is 
not generally known, an empirical equation is developed that relates 
the outlet cooling-gas temperature to known temperatures, flow con- 
ditions, and geometry. The empirical equation is 


~k^c c p , c 
ZU 



_ 0.29 b 1,2 
G g — 


T, 


6,2 


-0.38 


T, 


c,l 


-0.23 


( 20 ) 


where 

b cooling-passage height, (ft) 

c_. « specific heat at constant pressure of cooling gas, 
P '° (Btu/(Ib)(^R)) 
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G c 

G S 

l 

T c,l 

T e,2 

u 


mass Telocity of cooling gas, (lb/(sec) (sq ft)) 
mass velocity of combustion gas, (lb/(sec)(sq ft)) 
length, (ft) 

inlet cooling-gas temperature, (°E) 
outlet cooling-gas temperature, (°R) 
effectiTe OTer-all heat- transfer coefficient 


The factor K in equation (20) is shown to be a function of the 
ratio of heat-transfer surface on the cooling-gas side S c to the 
heat-transfer surface on the combustion-gas side S g . For longi- 
tudinal fins of l/l6 inch thick Inconel spaced l/Z inch apart and 
having heights of l/4 to l/2 inch. 


K - 6100 



( 21 ) 


which reduces to 50,000 for an unfinned annular passage where 
S 0 /S g equals unity. 

A combination of equation (20) with the characteristic equa- 
tion (18) 


r g,2 - T c,2 _ 


_ T g,l “ T c,l 




ZIT 


g,l \ bG c G P,c/J 


IU 


bG~c 


c°p,c 


T 


8,2 “ T g,l 


IU 


(18) 


tG c c p,< 


giTes the unknown outlet cooling-gas temperature in terms of known 
temperatures, mass flows, and geometry. 

Temperatures calculated by this analysis and the selected 
assumptions show good agreement with corresponding temperatures 
measured on seweral experimental tail-pipe burners. Calculated 
cooling-gas temperatures at a station 0,7 of the distance from the 
flame holder to the nozzle outlet are within ±15° of the corres- 
ponding measured temperatures; and calculated inside-wall temper- 
atures are within approximately ±25° of the corresponding measured 
temperatures. The difference between the measured and calculated 
temperatures of the inside wall ranges from about 0.06 to 0.16 of 
the difference between the measured temperatures of the inside wall 
and the cooling gas. 
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The agreement shown applies only to combustion chambers having 
a combination of well-proportioned flame holders and an optimum 
fuel-injection system, which results in a complete filling of the 
burner cross section with flame and a high combustion efficiency. 

The analysis predicts temperatures that are too high for con- 
figurations in which the fuel is concentrated near the center of 
the combustion chamber, or that have no flame holder, or where the 
radial clearance between the flame holder and the inside wall is 
large enough to provide a relatively cool layer of gas between the 
frames and the wall. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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APPENDIX A 


SAMPLE CALCULATION 

The Inside-wall temperature Is calculated fear a typical set 
of conditions obtained on an experimental tail-pipe burner. Badi- 
ation from the combustion gas is assumed to be twice the nonlumi- 
nous radiation and the outside-wall heat-transfer coefficient is 
estimated from turbulent flow over a flat plate 20 feet long. 

Given conditions. - The following quantities are assumed to 
be known: 


Burner length, 2, in. 48 

Outside diameter of Inside wall, D^, in 31 

Inside diameter of outside wall, D 0 , in. 32 

Mass flow of cooling gas, W c , lb/seo 3.14 

Mass flow of combustion gas, Wg, Ib/sec 35.78 

Tunnel airspeed, Vq, ft /sec 140 

Tunnel air temperature, Tq, 504 

Tunnel static pressure, p^, lb/sq ft absolute 781 

Combustion-gas temperature Tg at nozzle outlet, ... 3363 

Static pressure of combustion gas at 48-inoh station, 

lb/sq. ft absolute 1240 


The following temperatures were measured on the 48- inch section 
of the combustion chamber: 


Measured 

tempera- 

Temperature 

(°R) 

ture 

Inlet station 

Outlet station 

to 

1488 

1538 

?c 

1711 

1712 

H 

1744 

2130 

T S 

1704 

a 3114 


a Taken as inlet combustion-gas temper- 
ature plus 0.85 of combustion-gas 
temperature rise to nozzle outlet. 


Hydraulic diameters and passage height. - The hydraulic diam- 
eter of the annular cooling passage Is obtained from equation (4) 



NACA EM E9L09 


25 



= 0.08333 ft 


and the hydraulic diameter of the combustion-gas passage equals the 
geometrical diameter Dj_ 

D g - Di = || = 2.583 ft 
The cooling-passage height is 

b = D c - 0.04167 ft 

Floy areas and mass velocities, - The flow areas of the cool- 
ing passage and the combustion chamber are, respectively, 

it (32 2 - 31 2 ) 

A c = J A 144 L ~ 0*3435 sq ft 

Ag = j 2.583 2 * 5.24 sq ft 
from which the respective mass velocities are 


G c = o T ^r fs = 9.14 lb/(sec)(sq ft) 


G g = - ~y| = 6.83 lb/(sec) (sq ft) 

Eadiation between walla. - The inside and outside walls of the 
cooling passage were Inconel, which had an emissivity e v of about 
0.7 for the surface conditions and temperatures encountered (ref- 
erence 8). The corresponding value of e* v is 0.85 (equation (6)). 

From equation (8), 


F 


A,E 



0.543 


Eadiation from combustion gas . - The nonluminous radiation 
from the combustion gas is calculated frcta data in reference 8 by 
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assuming complete combustion of a stoichiometric mixture of G ±2^-26* 
for vhich the ratio of partial pressures of the carbon dioxide and 
water vapor to the static pressure of the combustion products are, 
respectively, 

p c° 2 

= 0.1588 

P g 


and 


v E20 

P g 


= 0.1276 


The beam length L is taken as 0.8 of the combustion-chamber 
diameter (reference 8, p. 69). The product 


P S L 


0.8 X 1240 X 2.583 
2116 


= 1.21 foot atmosphere 


and 


Pqq^L => 1.21 X 0.1388 => 0.168 

P H Q L = 1.21 X 0.1276 = 0.154 
2 


■ 00 . 


L + 




,L = 0.322 


The partial pressures of the carbon dioxide and water vapor in 
atmospheres are 

1240 

P C0 = ’Ini X 0 - 1588 “ °* 081 
2 

%,o - MS x °- 1216 - °-° 75 


and 


V 


p h 2 o + p co 2 


= 0.48 
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As an initial estimate, let 

t ± 2 = 2150° E = 1690° F and 

from the given conditions, T g 2 

= 3114° E = 2654° F for which 

6 co 2 = 0,057 . - 

“COg = °‘° 89 

e H20 = ®»039 

“HgO = °* 067 

G 1 = 1.05 

c i £ h 2 o = 0,041 

°i cc h 2 o = 0,070 

sa 0 « 005 

Aa = 0.004 

€g = 0 • 093 

a g = 0.155 


Use of empirical equation. - From equation (20), 
0.84 

n 0.29 b^ „ 

"Sj2 


■u n _ /n \U.o4 -j 

® PjC cq QOOi ®. 1 p 0.29 b _ 

KT -50,000^-| G g — 


-0.38 -0.23 

■ L c,l 


( 20 ) 


= 50,000(|^j 0 ’ 84 x6.83 0 * 29 X— X3114 -0 * 38 X1711 -0,23 
= 5.21 


and from the zero curve of 


T S,l- T o,l 

T S,2 _T g,l 


in figure 5, 


T g.?.g. ~ T .?-7. 2 . 

T S,2 " T g,l 


0.91 


Therefore, 


T 0#2 = 3114 - 0.91(3114 - 1704) = 1830° E 


Correction of cooling-gas temperature for heat losses. - The 
outside-wall heat-transfer coefficient is estimated from equa- 
tion (24): 


h c = 0.64 


% 


0.3 


(VqP) 


0.8 


■m 


. 0.2 


(24) 
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v here 


+ T, 


% 


■m 


1 ^ 1 ± 88 . , . + - 1 - 5 . 5 - 8 + 504^ = 1008° R 


and 


781 


53.3 X 1008 


= 0.0145 lb/cu ft 


The value of l, taken as the length of engine plus tail-pipe 
"burner, is 20 feet. 

h Q = 0.64 x 1008 0 * 3 X X ■ 4 * 95 Btu/(hr)(sq ft)(°R) 


From equation (22) 


a 0.173 XlO” 8 X 0.7 


S, 


f 1488 + 1538V 
\ 2 


- 504* 


1538 


4.95 (^ 1488 - 504 

= 11,260 Btu/(hr)(sq ft) 
and from equation (25) 


AT r 


11,260 X * X |§ X 4 


119° R 


° 3600 X 3.14 X 0.28 

The corrected cooling-gas temperature is 

T c#2 - 1830 - 119 = 1711° R 

Seat-transfer coefficient. - Equation (2) gives the heat- 
transfer coefficient for the cooling gas as 

17 V 1 0.3 v Q "I yf 0.8 

h 0 - 0.378 X “= = 34.1 Btu/(hr)(sq ft)(°R) 

0.0833°' 


and for the combustion gas as 
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hg = 0,378 X “ 16.27 Btu/(hr)(aq ft)(°R) 

2.583 

Detemlnatlon of inside-wall temperature , - prom equations (16) and. (17), 

5114 4 X 0 , 175 XlO" 8 X 0 . 85 X 2 X 0 .093 + 3 114 X 16 . 27 + 1711 4 X 0 . 173 X 10~ 6 X 0 . 545 + 17 11 X 54 . 1 

16,27 + 34.1 


0,175 X10“ 6 (0.545 + 0,85 X 2 X 0.155) 2836 5 = Q ^ 
T * 16,27 + 34.1 ~ * 


for which figure 4 gives * 0.774 

The calculated inside-wall temperature is 


t 1?2 » 2836 X 0.774 * 2195° E 


which Is close enough to the Initial estimate of wall temperature that a reestimation is 
unnecessary. 


NACA EM K9L09 


30 


NACA EM E9L09 


APPENDIX B 

DERIVATION OF LONGITUDINAL DISTEIBOTION OF 
COOLING-GAS IEMPERATOEE 

The one-dimensional distribution of cooling-gas temperature 
will be derived for the case of constant over-all heat-transfer 
coefficient U, a linear combustion-gas temperature distribution, 
and no heat lost to the outside. With reference to figure 1, the 
general case when T g ^ and Tc,l are different is considered. 

It is then assumed that at some station upstream (where x = 0) 

Tg «= T 0 «a T c . The linear distribution in combustion-gas tempera- 
ture is expressed by 



(Bl) 


The heat balance across a small segment dx of the Inside wall 
gives 


w c°p,c^c - U(T g ' *c> ^i*^ 
Using the simplifying notation that 


B 


UnDi 


^o°p,o 

U 


(B2) 


(B3) 


equations (Bl), (B2), and (B3) combine to give 

P 


&(T g - T 0 ) (T g - T c ) 


dx 


B 


which has the solutions 


and 


(Tg,l " Tc,l) = BP \1 - e B , 

/ V 


(*g,2 - T Cf2 ) -BP \l - 


B 


(B4) 


(B5) 

(B6) 
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By rearrangement and subtraction of equation (B5) from equation (B6) 


I . J - % - los, (i - lcSe (x . 


However, 


(B7) 


_1_ _ l /_l\ Z 1 _ 2U 1 

BP = B \j3zJ “ B (T g , 2 - T 8#1 ) = bG 0 c p ^ 0 {T g)2 - T g#1 ) 

and equation (B7) can be rewritten as 


T g,2 ~ T o, 2 _ 
T S,2 " T g,l 


1 - 


T - T 

i - SA 

( zu \ 

rn . rn 

- 8,2 ^S,l 

\b G c c p,c/ 


zu 


e bCi-c c p,c 


ZU 


(B8) 


bG c Cp^ c 


which relates the outlet temperature of the cooling gas to U and 
the known variables. 
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(a) Notation for combustion chamber and cooling passage. 



(b) Notation for longitudinal fin* Af = Jrl> 
C = 2l'y s c = (s + 2 lf)l. 


Figure 1. - Schematic diagram of combustion chamber and 

longitudinal fin. 
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Figure 3. - Effect of assumed distribution of combustion-gas tem^ 
perature on temperature of inside wall and cooling gas* 
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Figure 6. - Continued. Determination of exponents in empirical 

equation (19). 
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(d) Determination of m. 


Figure 6. - Continued. Determination of exponents in empirical 

equation (19). 
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Burner length, l, ft 


(e) Determination of u. 


• - Continued. Determination of exponents in empirical 
equation (19). 
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Figure 6. - Concluded. Determination of exponents in empirical 

equation (19). 
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Figure 7. - Determination of mean value of K for annular cooling passage, {See table I foj* 

symbols. ) 
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1>1 ^ I If 7f 8 P 4L T c,l T g»l T g»2 

(ft) (ft) (ft) (ft) (ft) (ft) (lb/ (aoc|) (lb/aq (°R) (°R) (°R) 

(aq ft)) ft aba* ) 


< 1.667 0.02083 5 10 1458 600 500 3600 

V 1*667 . 04167 5 10 1468 500 500 3500 

A 1.667 . 02083 5 0.02083 0.0052 0.04167 10 : 1458 500 500 3500 

q 1.667 . 04167 5 . 04167 . 0052 . 04167 10 1458 500 500 3500 
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Figure 9. - Experimental tail-pipe burner. Flame-holder blocked area, 240.5 square Inches ; 

percentage blocked area, 29.9. 
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Figure 10. - Variation of inside-wall and cooling-gas temperatures 
with oorabust ion-gas temperature at 48-inch station of tail-pipe 
burner shown in figure 9. 
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Figure 11. - Comparison of calculated and measured inside-wall and 
cooling-gas temperatures at 48-inch station of tail-pipe burner 
shown in figure 9. 
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